Tracer dispersion in the turbulent convective layer. 
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Experimental results for passive tracer dispersion in the turbulent surface layer under convective 
conditions are presented. In this case, the dispersion of tracer particles is determined by the inter- 
play of two mechanisms: buoyancy and advection. In the atmospheric surface layer under stable 
stratification the buoyancy mechanism dominates when the distance from the ground is greater 
than the Monin-Obukhov length, resulting in a different exponent in the scaling law of relative 
separation of lagrangian particles (deviation from the celebrated Richardson's law). This conclusion 
is supported by our extensive atmospheric observations. Exit-time statistics are derived from our 
experimental dataset, which demonstrates a significant difference between tracer dispersion in the 
convective and neutrally stratified surface layers. 

PACS numbers: 47.27.nb, 47.27.eb, 64.60.al, 05.45.Tp, 92.60.Fm, 92.60.Mt 



I. INTRODUCTION 

Understanding transport properties of turbulent con- 
vective flow (e.g. transport of particles, chemical species, 
temperature, etc.) is of significant importance for a num- 
ber of fields of science and technology covering many 
physical applications including: environmental pollutant 
dispersion, heat transfer, combustion, reacting flows, me- 
teorology and oceanology. Due to the buoyancy effect 
and the associated anisotropy, convective turbulence pro- 
duces much more complex phenomenology than the clas- 
sical Kolmogorov model of isotropic turbulence and this 
imposes new challenges for its analytical treatment 0- 
3]. In spite of this complexity there has been remark- 
able progress in understanding this phenomena based on 
the application of modern methods of theoretical physics 
(see for comprehensive reviews). Examples include 
the derivation of Kolmogorov and Bolgiano spectra in 
thermal convection, scaling of statistical moments of con- 
centration, multi-fractal structure of tracer statistics and 
some others (see [3, 17rll0j and references therein) . 

An important step for further refinement of the devel- 
oped theoretical framework is experimental validation of 
some of these theoretical predictions, which was motiva- 
tion for the reported study. In this paper we present our 
experimental results on the dispersion of passive tracers 
in the turbulent convective layer. The wall-bounded tur- 
bulent flow under the effect of changing stratification pro- 
vides flexible settings to study the relative contribution 
of the buoyant and kinetic energy fluxes on scalar trans- 
port since this contribution can be easily controlled by 
varying the distance to the underlying surface (and pass- 
ing the threshold of the Monin-Obukhov length) |2j, [7], |9| . 
The results pre sented here are an extension of our pre- 
vious study [ll| conducted for the case of the neutrally 
stratified turbulent surface layer. 

Various mechanisms of turbulent dispersion can be 



characterized by the mean inter-particle distance 0- 
ilH: 



(R 2 (t)) ~ At" 



(1) 



where n is the power exponent of the particle separation 
law and is specific for a particular mechanism of turbulent 
dispersion; A is a scale-independent dimensional parame- 
ter. A particular value of parameter n can be associated 
with an energy injection mechanism of turbulence and 
can be used for universal characterization of associated 
dispersive properties. For instance, for tracer disper- 
sion by Kolmogorov (locally isotropic) turbulence, n = 3 
(Richardson's law) and for Bolgiano-Obukhov (buoyancy 
dominated) turbulence, n = 5 [l|, 0, H|. In the surface 
layer turbulence the dominant dispersion mechanism is 
associated with the strong gradient of dissipation energy 
leading to 'ballistic' separation, n = 2 [ill ]. 

Based on the above comments we can draw an impor- 
tant conclusion regarding the time evolution of parameter 
n in the turbulent surface layer. More specifically, we in- 
fer that the long-time trends of this evolution should ex- 
hibit a striking difference in the convective and neutrally 
stratified surface layer. Indeed, initially in both cases 
n = 3 (this is valid for the short-times which correspond 
to the short spatial scales, so turbulence can be consid- 
ered as isotropic). Then in the neutrally stratified tur- 
bulent surface layer n decreases, approaching the global 
ballistic asymptote n = 2 [ll[ . Contrarily, in the convec- 
tive layer n increases reaching its 'convective' asymptote 
n = 5. Since a relative increase (or decrease) of param- 
eter n is quite significant, this difference should eventu- 
ally manifest itself in experimental observations and the 
statistics of tracer particles. Observing this phenomenon 
was the rationale behind and one of the main goals of the 
current study. 

The exponent n can be estimated from time series data 
of particle concentrations, as there is a well-known con- 
nection between the power exponent of the particle sep- 
aration law and passive tracer statistics [1, [T^j : 
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S 2 {R)=2[F 2 {Q)-F 2 {R)]o,R 2 l n , 



(2) 



2 



where 

F 2 (R) = (C(x + R,i),C(x,i)> (3) 

is the pair-correlation function of our concentration lev- 
els, C. These expressions are used to estimate the value 
of n in our model. 



II. EXPERIMENTAL PROCEDURE 

The experimental procedure was fully described in our 
previous publications [U E3| • Consequentially, here we 
briefly present details that are specifically relevant for the 
reported set of new measurements. 

Measurements of particle concentration were per- 
formed at a given location within the Taylor hypothe- 
sis of frozen turbulence R = Ut, where the wind speed 
U = const. From ((T|), it follows that: 

S 2 (R) = 2[F 2 (0) - F 2 (R)] oc R 2 ' n oc {Ut) 2 ' n , (4) 

so single point time-measurements of concentration also 
show power-law behavior with the same characteristic 
exponents as the spatial measurements. From these ob- 
servations we estimated the value of parameter n in ((4]). 

A continuous (and relatively stable) influx of particles 
into the turbulent flow was supplied by anthropogenic 
activity in the surrounding urban areas (about 20 km in 
size). The observation tower was at a height of H = 12 
m and the tracer concentration was measured by means 
of scattered light intensity off sampled particles. Con- 
tinuous observation was undertaken over a period of 21 
days with an air sample reading taken every 5 seconds 
(30 liters per sample) . The instrumentation provided to- 
tal particle count as well as particle size distribution in 
the range of 1 — 10 uia. Local meteorological observa- 
tions (three conventional meteorological stations in the 
vicinity) were available to draw conclusions about local 
meteorological conditions of the surface layer (stability, 
profiles of wind speed and temperature) . Wind speed and 
temperature were also measured directly at the sampling 
point (H = 12 m). The Monin-Obukhov length L mo was 
estimated from these inputs [3 EE) ■ 

As usual in meteorological studies, we anticipated that 
the night time observations would correspond to the 
boundary layer under stable or neutral conditions. Con- 
trarily, for convective turbulence conditions to be extant, 
an intensive heat flux from the underlying surface has 
to occur and this, of course, should include the signif- 
icant effect of solar radiation. Therefore, for study of 
particle dispersion in the convective regime of the tur- 
bulent surface layer we focused our observations on hot, 
sunny days during the morning hours (i.e. the measured 
air temperature was always above 20°C). The data sets 
with negative values of L mo were selected with antici- 
pation that additional constraints |L mo | < H and low 
wind velocity ensuring convective dominated mechanism 
of particles dispersion at sampling location. 



III. RESULTS 

Undertaking correlation analysis of the time series of 
the concentration C(t) of single point measurements, we 
found that the minimum time span that provides reliable 
statistics of the process should be more than two hours. 
As the first step, we plotted S2(R) as a function of R = 
Ut on a log-log scale and evaluated the scaling exponent 
n from the expression Sa(iJ) oc {Ut) 2 / n predicted by ^ 
and ([JJ. Some examples of these plots are depicted in 
Figffl 

In this paper we report on the analysis of the data that 
corresponds to the convective conditions of the turbulent 
surface layer. Similar results observed for the neutrally 
stratified surface layer [ll[ are presented for comparison. 

The values of mean wind speed Ujj and Monin- 
Obukhov length L ma relevant to the data series depicted 
in Fig. [T]are listed in Table U where Uh = U at eleva- 
tion from the ground, H (12 m). The value of Monin- 
Obukhov length L mo was estimated from local meteoro- 
logical station data for each observation period [Til [l5j . 
As mentioned above, for the hight H > \L mo \ the tur- 
bulence in the surface layer is dominated by the convec- 
tive flux, so we expect that the asymptotic value n = 5 
should always emerge. Contrarily, for the measurements 
in the non-convective scenario we anticipated to recover 
the global ballistic regime (n = 2) that is valid for tracer 
dispersion in the neutrally stratified turbulent surface 
layer. The classical Richardson regime with n — 3 should 
be only a short-time asymptote in both cases. Particu- 
lar scenarios approaching 'global' asymptotes (n = 2 or 
ii = 5) may be rather convoluted due to the shear dis- 
persion mechanism, which can provide a profound contri- 
bution during intermediate-times [Tlj . However, at the 
long-time limit; parameter n will eventually tend to ei- 
ther asymptote if the aforementioned stability conditions 
return to the turbulence surface layer. 

TABLE I. Mean values for wind speeds Uh and Monin- 
Obukhov length Lmo as used in Fig. [1] 



Sample U H (m/s) L mo (m) 

1. 2.19 ±0.39 -9.36 ±3.71 

2. 1.20 ±0.11 -5.55 ±1.32 

3. 1.64 ± 0.23 -7.43 ±1.70 

4. 4.00 ±0.12 -9.71 ±1.87 

5. 4.39 ±0.32 -8.80 ±3.34 



The existence of two significantly distinct asymptotes 
for particle separation clearly emerges from our experi- 
mental data (see Fig. [TJ. More specifically, we observe 
that at the short-time intervals, t — > 0, the parame- 
ter n w 3 (i.e. the classical Richardson regime valid for 
isotropic turbulence), then it always reaches two signifi- 
cantly different limiting values (insets in Fig. [1} : n w 5 
(dispersion in convective turbulence) or n ss 2 (dispersion 
in neutrally stratified turbulence). The intermediate see- 
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FIG. 1. Structure functions of the tracer concentration 
S2(R) oc R a oc (UHt) a on log-log scales: top plot - convec- 
tive conditions; bottom plot - neutral conditions (as published 
in [llj for reference). Insets show convergence of the scaling 
exponent n of the mean inter particle displacement converges 
over time ({R 2 (t)} oc t n ). Each data set in the convective con- 
ditions plot corresponds to a maximum of four hours of obser- 
vations, with wind velocities Monin-Obukhov lengths outlined 
in Tabled] +> Sample 1; o, Sample 2; *, Sample 3; •, Sample 4 
and □, Sample 5. a = 2/3 is the Richardson regime, a ~ 1/5 
is the convective regime, a fa 2/5 is the velocity shear mech- 
anism and a = 1 is the ballistic regime. Global convergence 
to the convective regime (a — 1/5) is clearly visible. 



narios of approaching these asymptotes may exhibit non- 
monotonic behavior, but eventually either the convective 
or ballistic regime will dominate. 

In order to further reveal the difference between par- 
ticle transport in convective and neutrally stratified tur- 
bulent surface layers we employ the algorithm of exit- 
time statistics for the concentration time series. These 
are statistics of time intervals where a measured value of 
concentration exits through a set of thresholds, 5C 
By scanning the time series for a given threshold, one 



can recover a set of times Ti(SC), where the measured 
concentration reaches this threshold. This set can then 
be used to calculate the Inverse Structure Function (for 
details see jfj): 

E q (6C) ee (t*(«7)>. (5) 

where q is a statistical moment. 

A comprehensive analysis of the properties of the In- 
verse Structure Function can be fulfilled by applying the 
well-known multifractal approach 0, 8, 161. This results 
in the following scaling: 

Z q (SC) oc («7)x(«>, x (q) = xmn[(q + 3 - D(h))/h], (6) 

h 

where h is the index of singularity from the range 
[hmini hmax] such that t ~ (SC) h and D(h) is the fractal 
dimension of the set with a singularity h. It is worth not- 
ing that particular values h = 1/5, h — 1/3 and ft. = 1/2 
correspond to convective dispersion, the Richardson law 
and ballistic scaling respectively, as discussed above. 

The function (q + 3 — D(h))/h is assumed to reach 
its minimum at the upper boundary of the singularities' 
range 16], so we can set h — h max and write 

X(<?) = (q + 3 - D (hmax)) /h m ax- (7) 

Similar to particle dispersion in the neutrally stratified 
turbulent surface layer, this leads to the conclusion that 
x{q) is a linear function of q We observe that x(q) 
seems to follow the predicted linear trend for q > 1, but 
with a noticeably larger gradient than in the neutral sur- 
face layer case (1.88 and 1.29 for the stable and neutral 
conditions respectively) . 
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FIG. 2. Mean Inverse Structure Function for experimental 
data sets displayed in Fig. [T] (convective conditions: □; neutral 
conditions [ill]: A). Error bars correspond to ± mean stan- 
dard deviation. The thin dashed (convective) and thick dash- 
dotted (neutral) lines represent a linear best fit over 1 < q < 5 
predicted by © [l]|. 
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As explained in @, 0| , a non-linear response near the 
value g = 1 (if it exists) can be attributed to a contribu- 
tion of the slow (differentiable) components of turbulent 
motion. In our observations we found that this change is 
usually less profound for the turbulent convective surface 
layer, see Fig. [2] 

IV. CONCLUSIONS 

We presented experimental results for passive tracer 
dispersion in the turbulent convective layer when the dis- 
persion of particles was dominated by buoyancy fluxes. 



We found that our observations can be intrinsically ex- 
plained with a two-stage model of tracer dispersion. Dur- 
ing the first stage of separation turbulence can be consid- 
ered as quasi-isotropic and tracer particles obey the stan- 
dard Richardson model. During the next stage (the long- 
time limit) particle separation follows the 'fast' regime of 
convective dispersion (predicted by Boligiano-Obukhov 
scaling). This is markedly different to particle dispersion 
in the neutrally stratified turbulent surface layer where 
the standard Richardson regime is followed by the much 
'slower' (ballistic) rate of separation. These results may 
have implications for the study of transport phenomenon 
in various geophysical systems. 
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